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The exact solution of the spin-i Falicov-Kimball model on an infinite-coordination Bethe lattice is 
analyzed in the regime of "classical" intermediate valence. We find (i) either phase separation or 
a direct metal-insulator transition preclude intermediate valence over a large portion of the phase 
diagram and (ii) within the intermediate valence phase, only continuous transitions are found as 
functions of the localized f-electron energy or temperature. 
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Introduction The phenomenon of intermediate va- 
lence (IV) is seen in a number of rare-earth compounds 
(such as Ce, SmB 6 , SmS, YbBi 2 , EU3O4, etc.), where 
the average f-electron filling per ion becomes nonintegral 
jjj . These materials exhibit distinctive anomalies in ther- 
modynamic and transport properties accompanied with 
either continuous or discontinuous valence transitions as 
the temperature, pressure, or composition are varied. 

IV materials have a localized f-level that lies near the 
chemical potential for the conduction electrons and is 
broadened by either hybridization with the conduction 
band or by electron correlations. This large density of 
states acts as an electron resevoir and leads to both the 
IV phenomenon and the thermal and transport anoma- 
lies. Theoretical descriptions have focused on two differ- 
ent approaches: (i) the Falicov-Kimball model (FKM) [^J, 
where the IV arises from an ensemble average of states 
with different integral valence ("classical" IV) and the 
transitions are either discontinuous or continuous, being 
driven by the strength of the Coulomb interaction be- 
tween the f-electrons and conduction electrons, or (ii) 
the periodic Anderson model (PAM) || , where IV arises 
from a quantum-mechanical mixture of states with dif- 
ferent f-occupancy ("quantum-mechanical" IV) and the 
transitions are continuous, being driven by the hybridiza- 
tion between the f-electrons and the conduction electrons. 
Most materials fall into category (ii) and it is believed 
that the experimentally observed discontinuous transi- 
tions arise from an additional coupling of the electrons to 
the lattice vibrations. A realistic theoretical description 
should include the physics of both of these approaches, 
but has proven to be cumbersome to carry out. Here we 
focus on case (i) to see why few materials can be found 
that fit into this scenario (likely materials 13,0] include 
EU3O4, EU3S4, and Sm 3 S4). We discovered that the 
"classical" IV phase is often precluded by either phase 
separation or a direct metal-insulator transition, which 
we believe helps explain why few "classical" IV systems 
can be found in nature. 

It is important to examine the experimental differences 
between scenario (i) and (ii). The most obvious differ- 
ence is in the response to a magnetic field. The FKM 



always possesses full f-moments, so it displays a Curie- 
like susceptibility proportional to 1/T. If one includes the 
additional (superexchange and RKKY) interactions be- 
tween the moments, then the uniform susceptibility will 
behave like 1/(T + T*), which will not diverge for posi- 
tive T* . Such a system will likely order in a spin-density 
wave at some characteristic temperature, though. The 
susceptibility never saturates or displays a maximum as 
T is lowered in the IV phase of the FKM (it does in the 
metallic phase when f-electrons are not present at low T 
and become thermally populated at higher T ||). The 
PAM, on the other hand, displays Curie-like behavior at 
high T but either saturates, or displays a low-T maxi- 
mum because of the Kondo effect, and the screening of 
the local moments by the conduction electrons. RKKY 
interactions are also present and can lead to magnetic or- 
der, which makes a sharp differentiation between the two 
models more difficult. Most real materials that do not 
have long-range magnetic order, display a susceptibility 
that either saturates or has a low-T maximum, indicat- 
ing that scenario (ii) applies. It is the purpose of this 
contribution to explain why scenario (i) is so difficult to 
attain. 

Model The spin-i FKM consists of localized f- 
electronic states and a delocalized conduction band. 
There is an on-site Coulomb interaction (U > 0) between 
the localized f-electron and the conduction electron. The 
model neglects the hybridization of the localized f-states 
with the conduction band, and the valence transitions 
occur only when the thermodynamic occupation of the 
different electronic states changes under the variations of 
external conditions (such as pressure, temperature, etc.). 
The Hamiltonian of the model is 

+ 4r**flfw - + , (1) 

where d\ G (di a ) is the creation (annihilation) operator 
for a conduction-band electron of spin a at site i, tij 
is the hopping matrix between lattice sites i and j, fj a 
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(fia) is the creation (annihilation) operator for a localized 
electron with its site energy Ef, and Uff is the on-site 
Coulomb repulsion between f-electrons. Uff is large in 
real materials, so we choose Uff — > oo, and restrict the 
number of f-electrons per site to nf < 1. We choose the 
total number of electrons to satisfy ritotai = nd + n/ = 1 
to examine the IV phenomenon where each ion donates 
one electron to the system. When the f-level lies below 
the bottom of the conduction band, the system is an 
insulator (rid = 0, rif = 1); when the f-level lies above 
the middle of the conduction band, it becomes a metal 
(rid = 1, nf — 0); IV phenomena can only occur when 
the f-level lies inside the bottom half of the conduction 
band. In our calculations, we adjust a chemical potential 
fi to satisfy the constraint ntotai = 1; /i is pinned to Ef 
in the noninteracting IV regime, and the average f-filling 
becomes nonintegral. 

Methodology The FKM can be solved in the infinite- 
coordination-limit, where the local approximation be- 
comes exact, and the momentum independent irreducible 
self-energy has a functional form which explicitly 

depends on nf, U, and the local Green's function G(ui). 
jj],[5|](| We examine this model on the Bethe lattice, where 
the density of states for the noninteracting system be- 
comes semicircular with the band width At* — At\fZ. 
[We take t* as our energy unit (i* = 1)-] In this case, 
there is a cubic equation for G(u>) that determines the 
interacting density of states A(ui) = — -lmG(ui) for any 
given nf Therefore, when n to tai = 1, we solve the 
problem by minimizing the free energy Q F[rif, ritotai = 
1] as a function of rif (0 < rif < 1): 

F[n f ,n total = 1] = 2 J dee /(e) A(e) + E f n f 
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+ 2Tjde {/(e)ln/(e) + [l-/(e)]m[l-/(e)]}A(e) 

+ T[n / lnn / + (l-n / )ln(l-n / )-n / ln2] , (2) 

where /(e) is the Fermi distribution function. 

First, we construct the ground-state phase diagram as 
a function of Ef and U (see Fig. |l|). For large enough 
U 0, there are only two phases: a metal when Ef > 
— ■J- and an insulator when Ef < — In this limit, 
the system becomes effectively noninteracting, avoiding 
the energetically unfavorable double-occupancy of both 
d and f electrons on the same site, and does not display 
IV. For U small enough, there is a range of values of 
Ef, lying in the lower half of the conduction band, where 
the chemical potential is pinned at Ef, and the average 
f-filling is nonintegral. As the system changes from a 
metal to a homogeneous IV phase, the value of the rif, 
at which the ground-state energy F gs [rif, ritotai = 1] has 
its minimum, increases continuously from 0. Hence, the 
boundary between the metal and the homogeneous IV 
phases may be obtained from the following condition: 
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(3) 



rif=0,n tota i=l 



which determines when the metallic phase is no longer 
a local minimum of the free energy. The analytic form 
of F gs [rif —* 0, ntotai = 1] is found from a low-density 
expansion || f°r ^( w ) an d G(uj). After some tedious 
algebra, the boundary equation resulting from Eq. (|^) 
becomes 
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Direct numerical calculations, minimizing the free energy 
in Eq. (g) shows good agreement with Eq. (||). This 
result is valid for U < 1.84177; larger values of U have 
the direct transition from the metal to the insulator at 
Ef = — 8/37T, which occurs before the metal becomes 
locally unstable. 

A similar analysis cannot be performed when nf — > 
1, because there is a first-order transition between the 
insulator and a phase-separated state. In order to 
show this phase separation, we minimize the free energy 
F[nf, ntotai] with respect to nf for fixed ntotai and then 
determine the free-energy curve as a function of ntotai- 
A Maxwell construction is finally performed to deter- 
mine the convex hull of F and see whether or not the 
unit-density case is phase separated. In equations, we 
compare 



F avg = aF[nf otal ] + (1 - a)F[nf otal ] 
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1], where the superscript A indicates 
> 1. 



< 1 and B indicates nf otal 
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Valence Transitions We begin at the U = limit. In 
this case, IV phases can be found whenever — 2 < Ef < 
and the chemical potential for the conduction electrons 
is pinned at Ef yielding n& conduction electrons (0 < 
nd < !)■ The remaining electrons 1 — n& are f-electrons, 
and the average filling per ion will be noninteger. How- 
ever, because all of the f-electrons share the same energy, 
the ground-state energy of this configuration is degener- 
ate with any phase-separated mixture of states with dif- 
ferent f-electron fillings (such as the integer-valent states 
nf = 0, and nf otal = n d and nf = 1, and ng tal = n d + l) 
because the f-electron energy is linear in the f-electron 
filling. In order to determine what situation is favor- 
able as U increases from zero, we need to expand the 
ground-state energy in a power series through second or- 
der in U, and determine whether the homogenous phase, 
or the phase-separated states are lower in energy. Such an 
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analysis is tedious, but shows that the maximal phase- 
separated state (where = and nj = 1) is stable 
whenever Ef < —1. The rest of the phase diagram is de- 
termined numerically in Fig. 1. There are four different 
stable phases as a function of Ef and U: (i) the insulat- 
ing phase, where all electrons are in the f-level, rid = 
and rif = 1; (ii) the metallic phase, where all the elec- 
trons are in the conduction band, = 1 and nf = 0; 
(iii) the phase-separated state, which is maximal when 
the A and B phases that the system separates into have 
rij = and nj = 1 respectively, and is intermediate- 
valent when at least one state has nonintegral f-filling; 
and (iv) the homogeneous intermediate valence state, 
where n/ is nonintegral. Notice how both phase sepa- 
ration and the direct metal-insulator transition preclude 
IV behavior over much of the phase diagram. The phase 
separation is likely to become inhomogeneous charge or- 
dering in a real material, because the long-range Coulomb 
interaction (ignored in the FKM) will not allow the sys- 
tem to separate into states that have excess charge over 
a large volume, and it will break up into microscopic do- 
mains of the different phases. 

We examine this behavior in more detail in Figs. 2 
and 3. Fig. 2 shows the Maxwell construction for the 
free energy in a case where phase separation occurs. No- 
tice how the free energy becomes concave near n to tai = 1 
which illustrates how the phase-separated state is stabi- 
lized. Fig. 3 is a vertical slice through the phase dia- 
gram at U = 0.9. It displays the characteristics of all 
different phases. When Ef < —1.265 the system is an 
insulator; for —1.265 < Ef < —0.995 it is a maximal 
phase separated state; for —0.995 < Ef < —0.845 it is 
IV phase separated, for —0.845 < Ef < —0.598 it is in 
the homogeneous IV phase; and for —0.598 < Ef it is a 
metal. There is no phase transition between the maxi- 
mal and IV phase-separated states, rather it is a smooth 
crossover, occurring approximately at the position of the 
circles in Fig. 1 (our criterion for the crossover is when 
nf < 0.999). We found that within the homogeneous 
iV phase, the valence change as a function of Ef (and of 
T) was always continuous, and exhibited no first-order 
(discontinuous) transitions. The only first-order transi- 
tions occur in the phase-separated states. Hence, the 
only way to have discontinuous IV transitions arises ei- 
ther from the coupling of the electrons to the lattice, or 
involves a phase-separation transition rather than a pure 
IV change. (We should emphasize here that we have not 
examined all other fillings of the FKM, where some dis- 
continuous transitions could occur.) These results differ 
from those recently seen in the one-dimensional model 
, where discontinuous IV transitions are also possible 
(although that calculation has difficulty differentiating 
from a discontinuous transition and phase separation). 

These results summarize the behavior at T — 0. At fi- 
nite temperatures, all of these phases survive, and phase 
transitions occur between them as a function of temper- 



ature. This is, in fact, the best place to look for these 
first-order transitions experimentally — by measuring the 
specific heat as a function of T and looking for the large 
spikes near the first-order transitions. Cy is determined 
by numerically differentiating the entropy (including av- 
eraging for the phase-separated states). Typical results 
are summarized in Fig. 4, for all different ground-state 
phases. All of our theoretical results are displayed in the 
experimental results of EU3S4 [Q: there is a first-order 
transition at 160K from a homogeneous to an inhomoge- 
neous IV phase (with a sharp peak in the specific heat 
accompanied by a structural transition) followed by a fer- 
romagnetic transition below 3.7 K (which is expected in 
any classical IV system which has local moments). 

Conclusions We have examined the phenomenon of 
intermediate valence in the spin-one- half Falicov-Kimball 
model. This is a model that can only display "classical" 
IV phenomena, as the microscopic occupation of the f- 
electrons is always exactly zero or one, but IV can oc- 
cur from ensemble averaging, where the average f-filling 
becomes nonintegral. Such systems should display mag- 
netic susceptibilities that are Curie-like (or Curie- Weiss- 
like if additional magnetic couplings are added in) which 
is not what is seen in most IV compounds. Instead most 
experiments show a susceptibility that either saturates, 
or has a maximum as T is lowered, which can best be de- 
scribed by the Kondo effect, and models that include the 
hybridization between f-electrons and conduction elec- 
trons. We discovered a fundamental reason why "clas- 
sical" IV materials are more difficult to find experimen- 
tally: either phase separation or a direct metal-insulator 
transition preclude the IV state over a wide range of 
parameter space. We believe this result helps explain 
why nearly all observed IV materials require hybridiza- 
tion to describe them. We also found, that the FKM 
does not support discontinuous IV transitions. All va- 
lence changes, within the IV phase, are continuous as 
functions of Ef or T. The only ways to get discontinu- 
ous transitions IV is to either couple the electrons to the 
lattice, or to have a transition to a phase-separated state 
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FIG. 2. Maxwell construction of the free energy (solid line) 
as a function of the total electron density, which shows that 
the system phase separates into an A-phase {nf otal < 1, circle) 
and a B-phase (nf otal > 1, square) rather than remaining 
homogeneous with n to tai = 1. The dotted line that connects 
A to B is the convex hull, and the vertical dot-dashed line is 
a guide to the eye for n to tai = 1- Here, U = 1.5, Ef — —0.85 
and T = 0. 
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FIG. 1. Phase diagram of the spin-i Falicov-Kimball 
model when ntotai = 1 an d T = 0. The filled-symbols and 
the solid lines indicate analytic results and the dotted lines 
are fit to the numerical results (open-symbols). The open 
circles denote the crossover from maximal phase separation 
(below the circles) to IV phase separation (above the circles). 




FIG. 3. Electron density for all four ground-state phases 
at U = 0.9 in the ground-state phase diagram. The solid 
symbols show ntotai an d the open symbols denote nj for the 
phase-separated cases. The dot-dashed line is the average 
electron density, and the dotted line is the average f-electron 
density. The system is an insulator for Ef < —1.265, phase 
separates when —1.265 < Ef < —0.845. It has IV when 
—0.995 < Ef < —0.598 and becomes a simple metal for 
Ef > —0.598. The dashed line between the maximal and 
IV phase-separated states marks the approximate location of 
the smooth crossover. 
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FIG. 4. Specific heat for the various different values of Ef 
when U = 0.9. The sharp jump as the temperature decreases 
indicates the first-order transition (the phase separation) at 
the corresponding critical temperature. The triangles mark 
phase-separated states. 
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